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ABSTRACT 

This is a proposal for an experiment to measure the relative 
stopping powers of Be, Al, Fe, Cu, Ag, W, Pb, and U by magnetically 
analyzing the degraded proton energies of the 742-Mev proton beam 
of the Berkeley cyclotron after passage through the absorbing elements. 
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RANGE-ENERGY MEASUREMENTS FOR 74Z-Mev PROTONS 

John E, Ward 

Lawrence Radiation Laboratory- 
University of California 
Berkeley* California 

May 4* 1959 



L INTRODUCTION 



If a charged particle passes through a substance with sufficient 
energy that the valence - shell effects--i, e. * atomic bindmg--can be 
neglected, the average rate of energy loss due to ionization only is 
expressed by the well-known Bethe equation, ^ 



dE 

dx 



a 4 2 
4-rre z 



NZ 



mv 



In - In(l-f5 2 )-|? 2 



( 1 ) 



where e and m refer to the electronic charge and mass, ez is the 
charge of the incident particle* NZ is the number of electrons per 
unit volume of the stopping element, and I is the "mean excitation 
potential" of the atoms in the stopping element. 

The development of this formula makes use of the Born approxi- 
mation, which requires the orbital velocity of the electrons to be small 
compared with the velocity of the incident particle, A proton with a 
kinetic energy greater than Z Mev satisfies this condition, ^ 

It can be seen that the stopping power depends only on the 
velocity and charge of the incident particle and not its mass* and the 
constant of primary importance as far as the stopping substance is 
concerned is its mean excitation potential. The mean excitation 
potential is defined as the product of all excitation levels in the sub- 
stance, K oo^> eac h weighted exponentially by its respective oscillator 
strength* 2 



I = 
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•^ik 



oo 
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where 




This formula would offer a means of calculation of I, but, except for 
the lightest of elements, a very impractical one. I is usually taken 
as an experimentally adjustable constant. 

Bloch has shown, using the statistical model of the atom, that 
in the Thomas -Fermi approximation one has 



I = KZ, 



(3) 



3 

where K is an adjustable constant. Using the more accurate Thomas - 

4 

Fermi-Dirac models Jensen found 




( 4 ) 



The index 0 refers to isolated atoms. This function is plotted in 
Fig. 1 as the solid line Q 



lo/Z, ev 
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Z 



MU- 17451 



Fig. 1. Experimental mean excitation potentials and the 

statistical atomic model. O: i/Z of bound atoms meas- 
ured with 340-Mev protons and evaluated relative to 
1^ = 165 ev; •: Iq/Z of isolated atoms. A: i/ Z meas- 
ured with 10-20-Mev protons and corrected for inner 
shells. The corrections apparently underestimate the 
stopping -power deficiencies of inner shells at low 
energies. 
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II. EXPERIMENTAL 

Since at high energies the corrections applied to the Bethe 
formula for valence-shell effects are negligible, the 742-Mev protons 
from the Berkeley cyclotron seem to be a natural choice for the incident 
particles. In addition, this energy range would serve as a check on the 
340-Mev Bakker and Segre data. ^ 

Figure 2 is a reproduction of a slide used by Werner Brandt at 
the Gatlinburg Meeting on Penetration (September, 1958), and compares 
the effective mean excitation potentials (without corrections), meas- 
ured by Burkig and MacKenzie at UCLA in 1957, ^ with the fluctuations 
of the electron densities in the periodic system. One can notice slight 
but significant trends in the I values, which are parallel to the electron 
density wherever the accuracy in I is high enough. The second max- 
imum should have been shown clearly by Mo (42) and Rh(45), but un- 
fortunately these were the poorest foils and hence gave the least accur- 
ate data in the series. Lead showed an "anomalously 11 high value and 
needs an independent check„ Thorium, on the other hand, confirmed 
the expected trend. 

On the basis of this diagram and on the advice of Dr. Brandt 

the following metals were chosen as target elements: Be(4), Al(13), 

Fe(26), Cu(29), Ag(47), W(74), Au(79), Pb(82), and U(92). It was 

decided to use enough absorber to degrade the protons approximately 

15% in momentum or from 742 Mev to 57 0 Mev in kinetic energy. In 

each case an equivalent amount of absorber was to be used so that the 

results could be compared on a relative basis with aluminum as the 

standard„ (See Table I) Adequate information exists on Cu, Al, Ag, 

Pb, and Be to determine the amount of absorber required. Specifically, 

7 

the tables and graphs m UCRL-2426 (rev. ) were used. To determine 
the equivalent thickness of Fe, W, Au, and U a numerical integration 
of Bethe 1 s formula was performed using values of I previously de- 
termined by Bakker and Segre (see Appendix A). Clearly these values 
are not precise. To provide for a means of interpolation in order to 
use the precise equivalent amounts of absorbers for the relative 
comparison, an additional length of aluminum absorber is to be added 
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Fig. 2. Mean excitation potential compared with electron 
density. (From Werner Brandt, Gatlinburg Meeting 
on Penetration, Sept. 1958. ) 
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to each absorber in order to provide an additional point. This absorber 
"shim, fr 9 g/cm^ of aluminum, is to be placed first on the "upstream 1 * 
and then "downstream” side of the various absorbers and the results 
averaged to give the effect of adding this additional absorber in the 



middle. 


The equivalent values of the aluminum shim 


are listed in 


Table II. 














Table I 








Equivalent absorber 


thicknesses 




Absorber 


Density 


Length 






(g/ cm 3 ) 


(g/ cm 2 ) 


(cm) 


(in. ) 


Be (4) 


1.85 


85 


45.94 


18.09 


Al(l 3) 


2.70 


89 


32.98 


12.98 


Fe(26 


7.85 


91 


11.59 


4.56 


Cu(29) 


8.93 


101 


11.30 


4.45 


Ag(47) 


10.49 


110 


10.49 


4.13 


W(74) 


18.60 


115 


6.20 


2.44 


Au(7 9) 


18.88 


118 


6.24 


2.46 


Pb(82) 


11.34 


133 


11.73 


4,62 


U(92) 


18.70 


125 


6.68 


2.63 


Table II 


Equivalent value of A1 shim 


Absorber 




EquiVo value a 


Thickness (g/cm^) 


Be(14) 




1.024 




9.216 


Al( 1 3) 




1.000 




9.000 


Fe (26) 




0.906 




8.154 


Cu(29) 




0.87 5 




7.875 


Ag(47) 




0.789 




7.101 


W (7 4) 




0.680 




6.120 


Au(7 9) 




0.676 




6.084 


Pb(82) 




0.660 




5.940 


U (92) 




0.630 




5.67 0 


a 

mass stopping power 
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The experimental setup is shown in Fig, 3, The full energy of 
the 184-inch cyclotron is used, 74Z Mev. The proton beam passes 
through a pair of jaws reducing it to l/l6 in, in the horizontal dimen- 
sion, The beam then passes through focusing and bending magnets 
and comes to a focal point 1 foot within the physics cave. This point 
is the physical center of the absorbers. The beam will then pass 
through a focusing magnet (three -element quadrupole), and a bending 
magnet in which the beam will be momentum-analyzed, and then fall 
upon x-ray film placed at the focal point of the quadrupole. 

The procedure will be to obtain an x-ray photograph of the full 
beam and then a photo of the beam with the absorber in place for each 
of the elements listed. In addition, photos will be taken with the 
aluminum shim as previously explained. 

The density of the developed x-ray film will be proportional to 
the number of protons striking it. By reading the film with a den- 
sitometer one can obtain a plot of the numbers of protons versus their 
positions on the film. In order to facilitate density measurements a 
duplicate of each of the exposures mentioned above will be made at 
precisely one-half the time for the original. Then, for example, the 
maximum intensity on the second exposure will be exactly one -half 
that of the maximum intensity on the first, and can be used as the 
reference for finding the points of half-maximum intensity on the first 
film. This procedure can then be extended ad infinitum throughout 
the intensity spectrum. The plot of numbers of protons versus film 
positions is expected to look like Fig, 4, In addition, by the use of 
wire orbits, a plot of energy versus position on film can be made, This 
plot should look something like Fig. 5. The two graphs can then be 
combined to form a plot of number of protons versus energy. An ex- 
ample of this type of plot is shown in Fig. 6. In general it has the same 
shape as Fig. 4. This plot (Fig. 6) will be the tool for describing the 
experimental results. 
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Fig. 3. Experimental setup, showing placement of magnets. 
Not to scale. 



Number of protons (orbitrory units ) 
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Position on film ( in ) 



MU- 17456 



Fig. 4. Expected general contour of plot of the number of 
protons, i. e. , the film density, versus position on an 
x-ray film after undergoing magnetic analysis. Two 
curves are shown: that of the full-energy proton beam 

of the 184-inch cyclotron, 742 Mev; and that of the 
same beam after approximately 15% momentum 
degradation. 



Energy ( Mev) 
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V 




Position of fi Im ( in.) 



MU- 17455 



Fig. 5. Expected general contour of the results of a wire- 
orbit analysis of an analyzing magnet. 
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Energy ( Mev ) 



MU- 17454 



Fig. 6. Expected results of the combination of data ob- 
tained from figures 4 and 5 showing the energy spectrum 
of the full-energy proton beam, 74Z Mev, and the same 
beam after approximately 15% momentum degradation. 
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III. DATA ANALYSIS 

For determination of the mean energy loss through each 
absorber the center of gravity of the two curves of Fig* 6 must be 
found. This can be done easily for the first curve, representing the 
full beam of protons undegraded. However, because of close collisions 
of protons and electrons the second curve has, theoretically, an 
infinite tail. Obtaining this infinite tail is not practicable experimentally, 
so that it must be supplied by a theoretical calculation. The observed 
ranges of individual particles from any monoenergetic source will show 
a substantially normal distribution about the mean range. Because the 
very hard collisions are few in number, the actual distribution is some- 
what asymmetric, with a long tail in the direction of the short ranges, 

g 

and a mean range which is slightly shorter than the modal range. For 
the approximation of this tail, the collisions suffered by the protons 
with the electrons of the absorbers will be divided into two classes, 
"near 11 and "far* " The so-called near collisions are those responsible 
for the occasional large energy losses and the far collisions are those 
contributing only to the mean energy loss, but introducing negligible 
energy-loss fluctuation. In other words, each proton will suffer 
approximately the same minimum energy loss and may in addition 
suffer a near collision (or collisions) resulting in a much greater energy 
loss. Each absorber will then be divided into a number of equal seg- 
ments in which the probability of a near collision is known, the energy 
loss will be computed for this one segment, and the resulting spectrum 
will be "folded" into itself as many times as there are segments. It 
is hoped that this spectrum will fit the experimental curve and provide 
the needed tail. This procedure seems valid in that the so-called near 
collisions are those in which the binding energies of the electrons are 
insignificant, and this is the assumption made originally by neglecting 
the corrections to Bethe' s formula. 

By this means of completing the graph for the degraded protons, 
the center of gravity of the curve can be found and the mean energy loss 
determined. This procedure would be done for each of the elements. 

By interpolation between the two curves, obtained by use of the 
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aluminum shim, for each of the elements the thickness of absorber 
equivalent to aluminum will be determined*, The relative mass stopping 

power of each element will be calculated by dividing the amount of 

Z Z 

aluminum (in g/cm ) by the amount of the element (in g/cm ) needed 

to produce the same mean energy loss. 

In addition, the stopping power per electron is, from Bethe' s 

formula, 



dE 

dX 

NZ 



. 4 Z 

4 ire z 



mv 



in 



Zmv 



in(l-p Z )-p 2 . 

_ 



The relative stopping powgr per electron with respect to aluminum is 



- ini + [in Zmv 2 - in(l-p 2 )-p 2 ] 

- inI Ai + [ inZmv 2 - in(l-p 2 )-p 2 ] 



This quantity can be easily calculated by multiplying the relative mass 
stopping power by 

A z Ai 

A Ai Z ' 

Thus the means for calculating I is apparent*. In this calculation the 

9 

value of I = 150 ev for aluminum is used. 

From the values of I calculated, a curve of i/Z versus Z will 
be plotted. This curve will in general look like Fig, 1, The value of 
the Bloch "'constant" can then be compared with previously obtained 
values. 

In order to determine the length of absorber traversed by the 
protons a correction for multiple scattering must be made. In the 
usual approximate correction, the observed range is expressed as a 
projection of the mean range: ^ 



R , 
obs 



y i. cos e. ~ y i,(i 4 e z ) 

t— i i Z_ i z i ' 

i i 



= R 

mean 



_1 

Z 
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where l is the distance between the ith and the (i-l)st small- 
1 — 

angle collision and 0^ is the direction with respect to the beam axis 
after the ith collision,, Therefore we have 




e s 2(x 

normal direction and x is the projected distance in the absorber. 
Also* 



is the mean square deviation in angle from the 



where 



< 



Lr ^ ( e s Z w) 



21 Mev i — r — 

o V t , 

pep 



where t is the absorber thickness measured in radiation lengths. 

12 

Using the value 





8 TT”Ne^Z 2 G(pv) 2 dx, 



Bichel 9 Mozley, and Aron found 



AR 

R 



= (R -R , )/R 

mean obs ' mean 



Gm^c^Z 

(M q c 2 ) 2 





0 - 1(3 2 

Rb7 z) 



dE 




Where G = 



2in 181 




B is the stopping number of the atom, 



B = Z 



Zmv 




and it 1 q is the rest mass of the electron 0 The ratio AR/R is a slowly 
varying function of the energy, amounting to several percent for high Z. 



-17- 



The corrections found by these experimenters appear in Table III. In 
each case it is clear that as Z increases the correction for multiple 
scattering increases, and as E increases this correction decreases. 



Table III 





Multiple - scattering 


corrections 


Element 


E 


Range obs a 




(Mev) 


Path length 


Be 


17.340 


0.9990 




9.588 


0.9988 


A1 


17.836 


0.9962 




14.971 


0.9960 




11.820 


0.9956 




6.150 


0.9943 


Cu 


17.893 


0.9900 




9.938 


0.9872 


Ag 


17.923 


0.9818 




10.022 


0.97 59 


Au 


17.549 


0.9626 




9.698 


0.947 5 



a 



correction for multiple scattering 
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